Estimating regional winter wheat yield by assimilation of time series of HJ-1 CCD NDVI into WOFOST-ACRM model with Ensemble Kalman Filter

Regional crop yield prediction is a significant component of national food security assessment and food policy making. The crop growth model based on field scale is limited when it is extrapolated to regional scale to estimate crop yield due to the uncertainty of the input parameters. The data assimilation method which combines crop growth model and remotely sensed data has been proven to be the most effective method in regional yield estimation. The methods based on cost function are powerless with crop dynamic growth simulation and state variable dynamic update. However, sequence assimilationmethodhas more advantages to overcome these problems, this paper presents amethod of assimilation of time series HJ-1 A/B Normalized Difference Vegetation Index (NDVI) into the coupled model (e.g. WOrld FOod STudies (WOFOST) crop growth model and A two layer Canopy ReflectanceModel (ACRM) radiative transfer mode) for winter wheat yield estimates using Ensemble Kalman Filter (EnKF) at the regional scale. The WOFOST model was selected as the crop growth model and calibrated and validated by the field measured data in order to accurately simulate the state variables and the growing process of winter wheat. The theoretically optimal time series LAI profile was obtained with the EnKF algorithm to reduce the errors which existed in both time series HJ-1 CCD NDVI and WOFOST–ACRM model. Finally, the winter wheat yield at the county level was estimated based on the optimizedWOFOSTmodel running on the wheat planting pixel. The experiment illustrates that in the potential mode, the EnKF algorithm has significantly improved the regional winterwheat yield estimates (R2 = 0.51, RMSE = 775kg/ha) over theWOFOST simulation without assimilation (R2 = 0.25, RMSE = 2168 kg/ha) at county level compared to the official statistical yield data. Meanwhile, in the water-limited mode the results showed a high correlation (R2 = 0.53, RMSE = 3005 kg/ha) with statistical data. In general, our results indicate that EnKF is a reliable optimization method for assimilating remotely sensed data into the crop growth model for predicting regional winter wheat yield. Crown Copyright© 2012 Published by Elsevier Ltd. All rights reserved.